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Microwave near-field quantum control of spin and motional degrees of freedom of 25 Mg + ions can 
be used to generate two- ion entanglement, as recently demonstrated in Ospelkaus et al. [Nature 
476, 181 (2011)]. Here, we describe additional details of the setup and calibration procedures 
for these experiments. We discuss the design and characteristics of the surface-electrode trap and 
the microwave system, and compare experimental measurements of the microwave near- fields with 
numerical simulations. Additionally, we present a method that utilizes oscillating magnetic-field 
gradients to detect micromotion induced by the ponderomotive radio-frequency potential in linear 
traps. Finally, we discuss the present limitations of microwave- driven two-ion entangling gates in 
our system. 



I. INTRODUCTION 

Scaling to large numbers of quantum bits (qubits) is 
important for the development of a future quantum in- 
formation processor. Several pathways toward this goal 
have been proposed based on using trapped atomic ions 
as qubits [IHB]. One particular approach envisions pro- 
cessing, transporting, and sympathetic cooling of ions in- 
side a large array of traps [2H3] . Recent experiments with 
various trap architectures have demonstrated the key in- 
gredients for scalable ion loading, trapping, processing, 
and transport [7HT7]. 

Current implementations of quantum information pro- 
cessors require a reduction of gate errors to be effi- 
ciently scalable [18], necessitating stringent qubit control. 
In most trapped-ion quantum information experiments, 
control is accomplished via laser-based techniques [4] , but 
other approaches are being investigated, e.g., those based 
on magnetic fields [T9H2T] . These methods reduce the re- 
quired laser overhead and may enable a higher level of 
control and integration than laser-based control. While 
the highest entangled-state fidelity is currently obtained 
with laser-based operations [22], magnetic-field based 
control may eventually provide sufficiently low errors for 
fault tolerant operations [23-27 . For single-qubit oper- 
ations performed with microwave magnetic near fields, 
errors per computational gate of 2 x 10 -5 and better 
have been reported [28j [29]. Although the theoretical 
threshold for fault tolerance depends strongly on the as- 
sumptions of the model under study, an operational error 
of 10 -4 or below may be a reasonable starting point in a 
future fault-tolerant architecture. Reaching this level of 
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precision for two-qubit operations requires improvements 
on the state-of-the-art errors by approximately two or- 
ders of magnitude. 

As with conventional computers, microfabrication may 
be advantageous in realizing a future quantum processor, 
such as those based on surface-electrode ion traps [30] 
(for a review see [31]). In addition to reducing the over- 
all volume of a processor for a given number of qubits, 
a reduction in scale could allow for tighter confinement 
with trap frequencies scaling approximately as 1/d 2 if all 
other parameters are held constant. Here d denotes a 
characteristic length scale of the trap, often character- 
ized by the smallest ion-electrode distance. Higher trap 
frequencies can enable shorter durations for transport, 
separation and recombination, and two-qubit operations 
mediated by ion motion. In the last two decades, typical 
trap length scales have been reduced from centimeters to 
tens of micrometers. This increases the oscillating mag- 
netic near-fields (scaling as originating from oscil- 
lating currents running in one or multiple trap electrodes. 
Fast and precise single-qubit rotations (carrier transi- 
tions) have been demonstrated in this way [28, 29} [32]. 
In addition, the corresponding field gradients (propor- 
tional to 1/d 2 ) can become large enough to couple the 
ions' spin with their motional state (sideband transi- 
tions) on experimentally accessible timescales [2TJ [32] . 
For trapped ions close to their motional ground states, 
the ratio of sideband-transition rates to carrier-transition 
rates driven by oscillating near fields and their corre- 
sponding field gradients is on the order of ao/d, where 
ao is the harmonic-oscillator ground-state extent [2T] . 
For typical experimental parameters, ao is on the order 
of 10 nm, and therefore ao/d ranges between 10 -3 and 
10 -4 for the smallest traps today. Further reduction in 
trap size would increase this ratio, enabling faster spin- 
motional entanglement. For the same reason, near-field 
gradients are typically not useful for larger (d > 100 jam) 
trap structures. 

The main practical limitation to size reduction has 
been "anomalous" motional heating, caused by electric- 
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field noise at the positions of the ions in excess of Johnson 
noise from the resistances of the electrodes and of their 
filters [33] [34]. Anomalous noise fields have been phe- 
nomenologically modeled as originating from potential 
fluctuations of surface patches or dipoles with typical di- 
mensions much smaller than <i, which yields a noise spec- 
tral density that scales as l/d A [33]. As d is decreased, 
the resulting heating grows faster than trap frequencies, 
oscillating near fields, and their corresponding gradients. 
Consequently, this form of heating appears to prevent 
scaling to even smaller trap dimensions. However, re- 
cent studies have demonstrated a substantial reduction 
of anomalous heating by cooling the electrodes [35H57] 
and by treatment of the electrode surfaces [38 , 39 . Such 
methods may enable the further reduction in trap size. 



In the experiments described here, we investigate 
building blocks for oscillating magnetic near-field quan- 
tum control of trapped ions [21]; a similar approach is 
pursued in [40]. To test the basic scheme, we integrated 
three microwave current- carrying electrodes into a room- 
temperature surface-electrode trap with d ~ 30 jam. Fre- 
quencies of the microwave currents are tuned to near res- 
onance with hyperfine transitions in 25 Mg + (nuclear spin 
5/2). The resulting magnetic field at the position of the 
ions determine the carrier coupling, while the magnetic- 
field gradient couples the hyperfine levels (spin states) to 
the motional degrees of freedom [21]. We previously re- 
ported on fast-carrier transitions and generation of two- 
ion entanglement [32 ; the entanglement was produced 
with a two-qubit gate operation [41-43 suitable for scal- 
able quantum information processing. In this report, we 
give a more detailed account of our experimental setup 
and the procedures for calibrating and optimizing the 
operation of the system. 
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In Sec. [IT] we discuss design, fabrication, and perfor- 
mance of the surface-electrode trap and describe the laser 
beams for loading, state preparation, and detection. Sec- 



tion III describes the microwave system that generates 
and shapes the signals used for coherent single-qubit and 
two-qubit operations. In Sec. |IV[ we give an experimen- 
tal sequence for calibrating and studying entangling-gate 
operations. The experimental procedure to adjust the 
microwave field at the ion's position for driving sideband 
transitions, a prerequisite for entangling gates [4TH43] . 
is described in Sec. IIV Al Section IIVBI details a tech- 
nique for rf micromotion detection by observing transi- 
tions in microwave magnetic-field gradients. In Sec. |IV C] 
we describe a procedure to optimize the orientation of the 
radial-motional modes relative to the microwave field gra- 
dient components and a static homogeneous quantization 
field. Finally, in Sec. [V] we discuss the factors currently 
limiting two-qubit gate fidelity [32] . 



FIG. 1. (Color online) (a) Cross-section view of the filter 
board and trap chip assembly, (b) Photograph of filter board 
and mounted surface-electrode trap, (c) Micrograph of the 
central region of the surface-electrode trap. The electrodes 
are labeled. The direction of the external quantization field 
| Bo | — 21.28 mT, located in the y-z plane, is shown. The 
25 Mg + ions are trapped at a distance d ~ 30 /im above the 
surface, at the position labeled trap center. 



II. SURFACE-ELECTRODE TRAP AND 
EXPERIMENTAL SETUP 

The surface-electrode trap is designed as an asymmet- 
ric "five-wire" trap [44] 05]. It consists of two radio- 
frequency (rf), six control, and three microwave elec- 
trodes; see Fig.[l]for a view of the trap assembly (FigJTk, 
b) and the central region of the trap electrodes (Fig. PTE) . 
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The trap chip is fabricated in the NIST Boulder clean- 
room. We choose aluminum nitride (A1N) as a sub- 
strate for its high thermal conductivity of approximately 
170 W/(m-K), which aids dissipation of heat generated 
by the microwave currents in the trap electrodes. Ad- 
hesion (Ti, 10 nm) and seed (Cu, 60 nm) layers are 
deposited on a 75 mm diameter A1N wafer by electron 
beam evaporation. A 14 /im thick layer of photoresist 
(Shipley SPR220-7 [46]) is spin-coated onto the wafer, 
exposed through a patterned chromium mask by con- 
tact lithography, and developed. The photoresist pat- 
tern defines 4.5 /im wide gaps between electrodes in the 
central region of the trap; the gap width is increased to 
25 /im at distances ~ 2 mm away from the trap center. 
Gold is then electroplated, in a 140 mL sulfite gold solu- 
tion (Transene TSG-250 [46]), up to a 10.5 /im thickness 
onto individual approximately 22.5 mm x 22.5 mm dies. 
The use of a small electroplating volume allows us to re- 
place the solution for each electroplating process, thus 
reducing the possibility of contamination. After remov- 
ing the photoresist with acetone, the seed and adhesion 
layers are wet etched, and a total of nine (approximately 
7.5 mm x 7.5 mm) trap chips are obtained by dicing. To 
provide good thermal contact, the chip is bonded with a 
thin layer of silver-filled vacuum compatible epoxy (EPO- 
TEK® H21D [46J to a solid copper support (Fig. [T^l) . 
This support structure also holds a printed circuit board 
for in- vacuum filtering of the control potentials. We use 
low-pass RC filters with R = lkQ (ANAREN, part no.: 
H2B15081001F2LO [46]) and C = 820 pF (NOVACAP, 
part no.: 0504N821J101P [46]). In addition, the printed 
circuit board is used for connecting both rf electrodes to a 
resonant quarter- wave step- up transformer (Q res — 350, 
when loaded with the trap) [47] . Connections between all 
trap chip electrodes and the filter board are made with 
gap-welded gold ribbons with a 50 /im x 500 /im cross 
section and length < 2 mm. Each microwave electrode is 
connected to a microstrip line on the filter board, which is 
soldered to a SMA jack on the input end. The other end 
of each microwave electrode is shorted to ground at the 
edge of the chip via gap- welded gold ribbons. In- vacuum 
coaxial cables connect the vacuum feedthroughs with the 
SMA connectors. All of the assembly steps are performed 
in a cleanroom environment to minimize the presence of 
contamination, e.g., dust particles on the surface of the 
trap chip, which can generate unwanted conduction paths 
and stray electric fields, possibly enhanced when exposed 
to UV laser beams. 

An rf peak voltage Vrf — 15 V to 60 V at ^rf — 
2tt x 71.6 MHz is applied to both rf electrodes, provid- 
ing the x-z radial confinement of the ions at a distance 
d ~ 30 /im above the surface. The rf frequency is stabi- 
lized to the center of the quarter-wave resonance of the 
resonator with a feedback loop to the frequency mod- 
ulation input port of the rf generator. The six control 
electrodes (CI to C6) are biased between —10 V and 
+10 V to provide confinement along the trap axis. The 
potentials are applied by six digital-to-analog converters 
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FIG. 2. (Color online) Calculated trap potential 0t rap (in 
units of meV) as a function of distance from the trap center 
(rf null and axial minimum), which is located 30 /im above 
the surface. Applied potentials are given in the text, (a) 
Potential in radial plane at y — 0. The dot marks the trap 
center. The vectors indicate the calculated motional mode 
orientation. Here «lf denotes the angle between the LF- 
radial mode vector and the x axis, (b) Potential in the axial 
direction at x = z — 0. The contribution from the axial 
rf pseudopotential (see text) is scaled by a factor of ten for 
visibility. 



(DACs; for details see Sec. Ill) 



Single-ion motional fre- 
quencies range between /axial = 0.5 MHz and 2.0 MHz 
in the axial direction, while two non-degenerate frequen- 
cies between / ra diai = 3 M Hz an d 12 MHz aligned along 
adjustable directions (Sec. IV C) in the radial x-z plane 
are used. 

We calculate the trapping potential 0trap in the 
radial and axial directions by means of the gap- 
less plane approximation [45 ; </>trap is generated 
by combining the rf pseudopotential and the con- 
trol potential. With V^f — 35 V and ^ci-C6 = 
{-0.801,0.641,-0.801,0.750,-0.384,0.750} V, Fig. ^ 
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shows the calculated ^>trap(^ 9 ^) i n the radial plane 
around the trap center at x = y = z = 0, while Fig. [2Jd 
illustrates 0trap(y) m axial direction. The overall trap 
depth, limited by a saddle point in the x-z plane (Fig.[2^i), 
is ~ 13 meV for this choice of potentials. We denote the 
two radial center-of-mass modes by their relative frequen- 
cies: low-frequency (LF) and high-frequency (HF) mode; 
their orientations are indicated by arrows in Fig. [2^. The 
orientation of these modes is determined by the control 
potentials while the rf pseudopotential is cylindrically- 
symmetric in a small volume around the trap center. 
Due to the tapered geometry of the rf electrodes, we 
find a weak axial rf pseudopotential, with a curvature 
corresponding to a motional frequency /axial — 20 kHz. 
The axial confinement, however, is dominated by the ap- 
plied control potentials and its frequency is simulated to 
be /axial — 0.51 MHz. Additionally, we notice the pres- 
ence of significant stray potentials, which must be taken 
into account in the simulations to reproduce the obser ved 
motional-mode directions and frequencies (Sec. IV C). 

We employ a set of "shim" fields, designed through 
simulations to produce orthogonal displacements of the 
ion(s) along x, y, and z while approximately preserv- 
ing the control-field curvatures, for micromotion nulling 
(Sec. IV B ) and to obtain oscillating magnetic-field maps 
). Due to the limited DAC resolution (Sec. Ill), 
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the ion position can be controlled to a precision of 
{Ax, Ay, Az} ~ {16,37,5} nm for typical trap frequen- 
cies / rad i a i ~ 7.0 MHz and / axia i ~ 1.5 MHz. 

A static magnetic field of magnitude |Bq| — 21.3 mT 
is aligned parallel to the trap surface at an angle of 
— 15° with respect to the z axis (Fig. and provides 
the ions' internal-state quantization axis. At this field 
strength, the ground state 2 Si/2 \F = 3, mp = 1) = |^) 
to 1^ = 2,777^ = 1) = |t) hyperfine transition frequency 
/o — 1.686 GHz is first-order insensitive to magnetic 
field changes, while the quadratic frequency deviation is 
~ 233 kHz/(mT) 2 ; we use the \\) and \\) states as the 
qubit. Here, F is the total angular momentum and mp is 
the projection of the angular momentum along the mag- 
netic field axis. Figure [3] illustrates the ground-state hy- 
perfine structure, highlighting the qubit states. We track 
| Bq | -field drifts every 20 to 30 minutes with a Ramsey 
experiment on the |3,3) to |2,2) transition (excited with 
microwave fields), which has a magnetic-field sensitivity 
of ~ 19.7 MHz/mT and adjust a computer-controlled cor- 
rection coil current accordingly. We observe day-to-day 
|Bo|-field variations of less than 0.02 mT. 

Two overlapping <j + -polarized laser beams for Doppler 
cooling on the 2 Si/2 |3, 3) -^ 2 P3/2 |4, 4) cycling transition 
propagate parallel to the magnetic field. A first beam 
(BDD), approximately —240 MHz detuned from this 
transition, performs initial Doppler cooling and optical 
pumping. A second beam (BD) is red detuned from this 
transition by approximately half the natural line width 
(r/2 ~ 2n x 20 MHz) and is used for final Doppler cool- 
ing and state preparation into the |3, 3) level of the 2 Si/2 
ground state [48 . State preparation and cooling takes 
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FIG. 3. (Color online) Ground state hyperfine structure of 
25 Mg + (nuclear spin 5/2) at |B | ~ 21.3 mT. The microwave 
driven hyperfine transitions used to prepare, detect, and con- 
trol the qubit states are depicted as arrows. The cooling and 
detection laser beam is indicated as a (blue) wavy arrow. (In- 
set) Illustration of the qubit states (bold lines) and their ra- 
dial motional states (gray). The "red" and "blue" sideband 
transitions as well as the qubit transition are indicated with 



1 ms; during the final 50 /is, only the BD beam is applied. 
For state detection, the BD beam is tuned to resonance 
for approximately 200 /is, discriminating the |3,3) state 
from the other hyperfine ground states through resonance 
fluorescence. The fluorescence photons are detected by a 
photon-multiplier tube (PMT) detector. In most exper- 
iments, the 1 3, 3) population is transferred via two mi- 
crowave preparation pulses to \\) to initialize the qubit. 
To detect the \\) state, two pulses are used to transfer 
the |^) population to the |3,3) (bright) state, while the 
|t) population is "shelved" into the |2, — 1) (dark) state 
(Fig. [3]). During the detection period, the PMT registers 
^ 11 counts for the bright state and ~ 0.2 counts for 
the dark state on average per ion. Stray light accounts 
for ~ 0.1 counts on average. The BDD and BD beams 
are both focused into the vacuum chamber with a 10- 
cm focal length lens, and their powers are adjusted to be 
~ 1 /iW and ~ 30 /iW, respectively; both have waists [49] 
of ~ 15 /im at the position of the ions. To create ions, 
a third continuous-wave (CW) laser beam (denoted PI: 
A ~ 285 nm, power ~ 3 mW, and beam waist ~ 15 jam) 
is directed antiparallel to the cooling beams and is used 
to photoionize neutral Mg atoms inside the trapping re- 
gion. An oven tube, holding enriched 25 Mg, is mounted 
approximately 2 cm away from the trap center. After al- 
lowing the oven to heat up for 30 s, the PI is applied for 
a few seconds until ion(s) are loaded. The oven current 
is adjusted to a loading rate of less than one ion per sec- 
ond, so the PI beam can be turned off manually after the 
desired number of ions are detected on a CCD camera. 
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FIG. 4. Schematic diagram of the microwave system used to 
control the microwave fields applied to the ion(s). It enables 
individual control of phases and amplitudes of the currents 
in all three microwave electrodes for two frequencies / s = 
fo =b /radial (BSB and RSB) around 1.69 GHz, within a range 
of about 200 MHz. It also enables a signal between 1.2 GHz 
and 2.3 GHz applied to MW2 for hyperfine state control and 
single-qubit operations. The gray components are computer 
controlled. 



At a background pressure of 2 x 10 -9 Pa, the storage 
period of a single ion is observed to be up to a few hours, 
with Doppler cooling applied periodically every one or 
two milliseconds. Typical trapping lifetimes for two ions 
are on the order of 30 min. 



III. MICROWAVE SYSTEM 

To implement microwave control of the ion(s) spin and 
motional states, we built a system assembled mostly from 
commercially available discrete microwave components 
(Fig. [4|. For hyperfine state manipulation and single- 
qubit operations, magnetic fields from microwave cur- 
rents in electrode MW2 drive transitions at frequencies 
between 1.2 GHz to 2.3 GHz, spanning all hyperfine tran- 
sitions in the 2 Si/ 2 manifold. As noted above, we use a 
series of two transfer pulses (carrier 7r-pulses) to transfer 
population from |3,3) to |3, 1) via |2,2), to initialize the 
qubit state. We have demonstrated 7r-times for qubit spin 
flips as short as ~ 20 ns [32]; however, here we adjust the 
power of the microwave signals to levels corresponding to 
7r-times on the order of a few microseconds, limited by 
the timing resolution of our data acquisition system (see 
below) . 

To drive "blue" (BSB) and "red" (RSB) sidebands of 
the qubit transition (Fig. [3]), we apply currents simulta- 
neously to MW1, MW2, and MW3 at f s = f + / rad iai or 
fs — fo — /radial- The amplitudes and phases of the three 
signals are adjusted to minimize the oscillating magnetic 
field at the position of the ion(s), while maximizing the 



six phase shifters and six attenuators. 



The microwave setup is integrated with a data acquisi- 
tion system, controlled by a field-programmable gate ar- 
ray (FPGA) with a timing resolution of 16 ns [50 j. The 
entire system comprises 16 digital (TTL) output chan- 
nels, 24 DAC channels, eight direct-digital synthesizer 
(DDS) modules, and two digital (TTL) inputs for regis- 
tering PMT counts. Each DDS can generate signals at 
frequencies up to 400 MHz with a frequency resolution of 
0.233 Hz and phase resolution of 0.022°. The DACs have 
an update rate of 500 kHz and can be set between — 10 V 
and +10 V with a resolution of 0.305 mV. We use six 
of these channels to apply potentials to the trap control 
electrodes and 14 channels to control microwave com- 
ponents. In addition, a fast, two-channel DAC is used 
to generate arbitrary waveforms with a 50 MHz update 
rate, a voltage range from —10 V to 10 V, and a reso- 
lution of 0.305 mV. This DAC is programmed via USB 
and triggered by the data acquisition FPGA. 



We use three DDS modules as sources for three fre- 
quency octupling modules, each consisting of three dou- 
bling stages in series with amplifiers and filters. Two 
of these generate the RSB and BSB signals, which 
pass through individual pulse shaping and control stages 
(phase and amplitude control). The pulse envelopes for 
the BSB and RSB signals are produced with a multiplier 
(Analog Devices ADL5391 [46]) used to shape pulse am- 
plitudes and controlled by one channel of the fast DAC. 
In addition, a voltage-controlled attenuator adjusts the 
overall power of a pulse. After this pulse-shaping stage 
the signals are split and directed to three individual 
voltage-controlled phase shifters (~ 12° /V) and attenua- 
tors (~ 1.6 dB/V), to control the power and phase of the 
signals delivered to the individual microwave electrodes. 
Before the RSB and BSB signals are combined onto the 
inputs of each microwave electrode, they are amplified 
by high-power amplifiers (Mini-Circuits ZHL-30W-252- 
S+ [46]). Following three high-power combiners, three 
coaxial cables connect the microwave setup with the vac- 
uum feedthroughs. The third octupling module generates 
frequencies to drive the hyperfine carrier transitions, and 
its signal is combined with one of the sideband lines at- 
tached to electrode MW2. Its power is adjusted with 
fixed attenuators to yield carrier 7r-pulse durations of a 
few microseconds. 



On a daily basis, we calibrate the frequencies and n- 
pulse durations for the transfer and single-qubit pulses. 
From day to day, we observe fractional variations of 
< 10 -3 for the pulse durations. Due to our quantiza- 
tion field calibrations and adjustments, the hyperfine- 
transition frequency values stay constant within a sta- 
tistical uncertainty of approximately 0.3 kHz. 
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rf quadrupole electric field 




FIG. 5. (Color online) Calculated trap rf quadrupole electric 
field in the radial plane. The quadrupole axes are rotated by 
a-RF — —20°. The arrows indicate the horizontal (red) and the 
vertical (gray) micromotion component. Only the horizontal 
component can be detected by the detection laser beam. The 
red and gray lines indicate the two shim directions used for 
the compensation of horizontal and vertical micromotion. 



IV. EXPERIMENTAL PROCEDURES AND 
CALIBRATION SEQUENCES 

We use a single ion to calibrate experiments and set 
parameters for the two-qubit gate. After an ion is loaded 
into the trap, micromotion along the projection of the de- 
tection laser beam in the radial plane (z direction) is min- 
imized by measuring the variation in laser fluorescence 
from the BD beam [51] as a function of the "horizontal" - 
shim field (Fig. [5]). Since the trap rf quadrupole is tilted 
in the radial plane by cvrf — —20° from the x axis 
(Fig. [5]) , the horizontal-shim field is chosen to push along 
a vector rotated by 2cvrf from the z axis. A push along 
this axis will result in a maximal variation of the horizon- 
tal micromotion amplitude. We verify this compensation 
field with a nulling scheme usin g micr owave magnetic- 
field gradients (described in Sec. IV B). This technique 
can also be used to minimize the micromotion in the ver- 
tical (y) direction. We find that the trapping fields are 
perturbed by stray fields, with corresponding potential 
Astray The exact origin of these fields is unknown, but 
they are consistent with a contact potential patch, and 
can be compensated for as described in Appendix [Aj 

Following micromotion nulling, we calibrate the mi- 
crowave magnetic fields used for dri ving th e radial side- 
bands of the qubit transition (Sec. IV A). From these 
experiments, we extract the orientation and the strength 
of the microwave magnetic- field gradient. As a last cali- 



bration step to set up the two-qubit gate experiment, we 
align the radial-mode vectors relative to the microwave 
magnetic field by rotating the modes with an additional 
control field (Sec.|lVC|). 



A. Controlling and Mapping the Microwave 
Near-Field 

To implement motional sideband transitions at fre- 
quencies / s , a gradient of the microwave magnetic field 
is necessary. However, any oscillating field at the ion(s) 
position will cause off-resonant carrier transitions and ac 
Zeeman shifts [21 , both of which will inhibit precise con- 
trol. Thus, the magnetic field amplitude should be sup- 
pressed as much as possible at the position of the ion. We 
realize this experimentally by adjusting the phases and 
amplitudes of the three microwave currents to minimize 
ac Zeeman shifts detected by the ion. The experimental 
sequence is shown in Fig. [6^l. 

The nulling procedure relies on the measurement of 
ac Zeeman shifts on superpositions of selected hyperfine 
states I a) and \b). In a coordinate system with unit vec- 
tors e 5 = B /|B |, e% = e x and = e 5 x e 5 , the 
magnetic-field components oscillating at f s can be writ- 
ten as: 



B 



MW 



1 / cos /3 — sin /3 

: _ e 2nif s t gin/3 cos/3 Q 



Bj_e n cose 
B±ie n sine 
Bn 



(1) 



where Bn and B± are the strengths of the oscillating 
magnetic-field components parallel and perpendicular to 
the quantization field Bo- Here (3 is the rotation angle of 
the perpendicular magnetic-field component polarization 
ellipse with the x axis and 7 is a phase. For e = 0, the 
oscillating magnetic field in the plane perpendicular to 
Bo is linearly polarized. For e = ±7r/4, the perpendicular 
field component has pure circular polarization. The Bn 
field couples levels with Ara^ = 0, whereas the B± field 
drives transitions with Amp = ±1. These fields induce 
an ac Zeeman shift Sf\cz = c^B\ + c\\B% on the \a) <o> 
I b) transition, where the coefficients may depend on f s 
and e. Note that we implicitly include the Bloch-Siegert 
shifts (energy shift resulting from the counter-rotating 
part of the field) in the calculation of the coefficients in 
£/acz- 

First consider \a) = |2,0) and \b) = |3,0). A B± com- 
ponent applied near fo is red detuned from the |2,0) — >• 
{|3, -1) , |3, 1)} and |3,0) H> {|2, -1) , |2, 1)} transitions, 
thereby increasing the energy of \a) = |2,0) and de- 
creasing the energy of \b) = |3,0). From the Clebsch- 
Gordan coefficients for these transitions, we determine 
c± ~ (0.357 - 0.092 sin 2e)Hz/(/iT) 2 at / s = f . A B\\ 
component near fo is red detuned from the |3, 0) —> |2, 0) 
transition and also increases the \a) <H> \b) transition fre- 
quency: Cm ~ 0.468 Hz/(/iT) 2 at / s = f . Both c± and 
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FIG. 6. Realization of optimal microwave- field configura- 
tions for driving motional sideband transitions, (a) Pulse 
sequence to detect the presence of residual oscillating mag- 
netic fields at the ion(s) position through phase accumula- 
tion from the corresponding ac Zeeman shift, while the ion 
is in a superposition state of \a) and \b). The spin echo tt- 
pulse makes the result insensitive to magnetic-field fluctua- 
tions slower than the time scale of an individual experiment 
(typically c± 500 /is) . We observe minimal phase pickup when 
the oscillating magnetic field is minimized at the position of 
the ion. (b) Extension to measure the ac Zeeman shift as a 
function of displacement 5r (see text), (c) Experimental data 
mapping out the ac Zeeman shift in the radial x-z plane by use 
of the superposition state of \a) = |2,0) and \b) = |3, 1). The 
fitted isolines are obtained from the quadrupole oscillating 
magnetic field component B\\ along Bo through its resulting 
ac Zeeman shift £/acz, which is proportional to B^. The an- 
gle ojmw is the rotation angle of the microwave quadrupole 
field, but the directions of minimum and maximum values of 
\5B\\/5r\ are along angles 2ojmw and 2ojmw + 90° (Fig. |8^l). 



cm are positive, of similar magnitude for all e and depend 
only weakly on / s around fo within the useful range of 
motional sideband frequencies, so this choice of \a) and 
| b) can be used to minimize both B± and B\\. 

Another useful choice is \a) = |2, 0) and \b) = 
1 3, 1). Here, the dominant ac Zeeman shift results from 
B\\ coupling to the |3, 1) —> |2, 1) transition: cy ~ 
-49.08/((/ 8 - / )/MHz)Hz/(/iT) 2 , and c ± ~ (0.199 + 
0.220 sin 2e) Hz/(//T) 2 at f s = f . This transition can be 
used to minimize B\\ even further due to its large cy value, 
by use of currents in electrodes MW2 and MW3 (those 
electrodes contribute dominantly to Bu due to their ge- 
ometry). We produce the data of Fig. pp using this tran- 



sition and an experimental sequence presented in Fig.[B]3. 
With Bmw = at position 5r = {fe, Sy : Sz} = {0, 0, 0} 
in the coordinate system of Fig. [lb [not the frame of 
Eq. ([!])] we expect to first order 

B M w =B' cos(27r f s t) 

cos(2aMw) sin(2aMw) \ / Sx 
fo/ 

sin(2aMw) — cos(2aMw) / \ Sz 

" (2) 

Here B' and c^mw characterize the strength and orien- 
tation of the quadrupole for small displacements S r < 
1.5 /im from the field null point. By fitting the ac 
Zeeman shift resulting from the projection Bu of Bmw 
onto B to the experimental data of Fig. [6p, we ex- 
tract B' = 35.3(4) T/m and a M w = -26.5(7)°. From 
this fit, we determine a small shift of the center of the 
magnetic quadrupole with respect to the trap center by 
{5x,5z} ~ {—0.11,0.23} /im, limited by the sensitivity 
of our nulling procedure. We detect phase accumulations 
over a time T p = 250 /is; we cannot measure phase shifts 
smaller than 7r/10 during that time. This limits the ac 
Zeeman shift sensitivity to ~ 0.2 kHz and leads to finite 
fields at the trap center. The experimental day-to-day 
variations for the nulled configuration are < 0.016 dB for 
the relative powers of MW2 and MW3 with respect to 
MW1, and < 0.12° for the relative phases of MW1 and 
MW3 with respect to MW2. 

Numerical simulations of the microwave currents in the 
three electrodes by use of ANSYS HFSS software [46] 
provide a value of B' sim 1A ~ 44 T/m for a current am- 
plitude of 1.0 A in conductor MW1, while the ampli- 
tudes and phases of currents in MW2 and MW3 are 
adjusted to give a nulled field configuration at the po- 
sition of the trap rf null. We find the following cur- 
rent amplitude ratios R and phase shifts A<£ in MW2 
and MW3 relative to the current in MW1 for the 
nulled configuration: {Rmw2, Rmws} — {0.82,0.96} and 
{A$mw2, A$mws} - {161.6°, -14.6°}. A comparison 
with the experimental current ratios is difficult because 
of slightly varying attenuations in the signal paths lead- 
ing to the three microwave electrodes. However, from 
simulations, we extract an angle q^w — —25°, close to 



the value fitted to the experimental data (also Sec. IV B ). 
A similar level of agreement between simulation and ex- 
periment is found in [40 . In our experiments, current 
densities are simulated to be as high as ~ 2 x 10 10 A/m 2 
near the edges of the conductor MW1, with a total cur- 
rent of 0.8 A in MW1, for the above stated experi- 
mental settings, producing a gradient of 35 T/m. We 
estimate the corresponding power dissipation in all elec- 
trodes due to all currents from the bulk resistance of 
gold 2.2 x 10 -8 Q-m) in the central region of the chip 
(± 200 /im along y) to be ~ 80 mW. In Fig. [7] the current 
density distribution in electrode MW1 and the induced- 
current densities in the neighboring electrodes RF1 and 
C2 are shown. The simulated magnetic and electric mi- 



8 



microwave current 
density 




microwave fields 



(a) 
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FIG. 7. (Color) Cross-section view of the simulated mi- 
crowave current density distribution in the trap electrodes 
when the magnetic field is nulled at the ion position, (only 
showing electrodes RF1, MW1, and C2). The densities shown 
here correspond to a gradient of B' sim c± 35 T/m. Currents 
in microwave electrodes induce currents in neighboring elec- 
trodes, here illustrated for RF1 and C2. 



crowave fields for the nulled configuration are plotted in 
Fig. [8] The influence of the microwave electric fields on 
the trapping potential is discussed in Appendix [Bl 



B. Microwave Approach for rf Micromotion 
Nulling 

In the presence of excess rf micromotion, the mi- 
crowave magnetic near-field gradient can also drive 
micromotion-induced sidebands on carrier transi- 
tions [51 j. This effect can be used to null micromotion 
along any direction in the radial x-z plane (refer to 
the frame of Fig. [TJd). Therefore, it may have an 
advantage over certain optical schemes that are sensitive 
to micromotion only along the direction of a laser beam 
(or along the direction of the effective /c-vector for 
stimulated- Raman transitions) [51]. For example, in a 
surface-electrode trap, the direction perpendicular to 
the surface can be difficult to probe with laser beams, 
because scattered light from the trap surface may 
interfere with detection and/or cause electrode charging. 

The micromotion can be written as (xmm&x +yMM£y + 
zmm^z) cos(27r/Rpt). In the presence of a microwave 
magnetic field gradient of frequency / s and in the coordi- 
nate frame of Fig.[]J, the ion in its rest frame experiences 
fields with frequencies / s ± /rf and amplitude Bi on , cf. 
Eq. (§: 



l / cos(2a M w) sin(2a M w) 
-B' 

* sin(2aMw) — cos(2aMw) 





f %M \ 




i 







(3) 



For f s = /hfs ± /rf, where / H fs is the frequency of a 
transition between states in the hyperfine manifold, the 
ion experiences an oscillating field at frequency /hfs- 
The resulting Rabi rate is proportional to |B ion | and 
therefore to tmm, the micromotion amplitude. To min- 
imize the excess micromotion, we move the ion in the 
radial plane (along the shim directions illustrated in 
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FIG. 8. (Color online) Illustration of the simulated microwave 
fields in the radial plane, (a) Magnetic quadrupole field as a 
function of displacement from the trap center, (b) Electric 
field. Due to the pseudopotential from the electric field the 
ion(s) are displaced from the trap rf null as indicated in gray. 
The pseudopotential from the electric field displaces the ion(s) 
from the trap rf null (see Appendix |b| . The estimated ~ 0.4 
nm displacement for typical trap parameters is shown in gray, 
scaled by a factor 10 3 for visibility. 



Fig. |5| to a position that minimizes the Rabi rate of 
the micromotion sidebands. The oscillating field B ion 
has projections B\\ and B± on Bo- Therefore, we can 
sense micromotion along any direction in the xz plane, 
since it is always detectable by either a Atuf = or 
Am^ = ±1 transition (or both). In the experiment, 
we use the |2,0) O |3, 1) and |3, 1) O |2, 1) transitions, 
with a magnetic-dipole matrix element ~ 0.414 ji^ and 
~ 1.001 /iB, respectively, where [1b denotes the Bohr 
magneton. 

To estimate the sensitivity of this nulling method, we 
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z(|jm) 

FIG. 9. Rabi rate on the red micromotion sideband of the 
|2,0) <B> 1 3, 1 ) transition as a function of ion position. The 
isolines result from a model fit to the data, from which 
B' — 35.1(4) T/m and ckmw = —31.1(2)° are extracted from 
Eq. 



assume B' = 35 T/m and a transition matrix element for 
/hfs of //b- To further simplify the estimate, the rotation 
angles of the microwave-magnetic quadrupole (ccmw — 
—30°, Fig. [8^i) and rf-electric quadrupole (q;rf ~ —20°, 
Fig. [5| are assumed to be oriented such that in the rest 
frame of the ion, we achieve a maximum oscillating field 
component that drives the transition (c^mw = <^rf/2). 
For the given misalignment, this approximation is good 
to within 15 % for the estimated sensitivity. Note that 
in the experiments both components of the microwave- 
field gradient can be used, and the misalignment does 
not limit the sensitivity of this method. The Rabi rate 
for the micromotion sideband-induced carrier transition 
is then given by [21]: 



MM 



2 BrMM 2h' 



(4) 



For our parameters, and assuming a 7r/10 rotation de- 
tection sensitivity on the micromotion sideband during 
a drive time T p = 2 ms, a micromotion amplitude of 
fMM =0.1 nm can be detected. 

In addition to detecting micromotion, we use this 
method to probe characteristics of the microwave- 
magnetic and rf-electric field by measuring the micro- 
motion sideband Rabi rate as a function of position. 
Figure [9] shows the Rabi rate on the red micromotion 
sideband of the |2,0) |3, 1) transition versus ion po- 
sition. Using these data together with the calculated 
shape of the rf quadrupole, we extract the parameters 
B' and olmsn of the microwave quadrupole field, finding 
B' = 35.1(4) T/m and a M w = -31.1(2)°. This map 
was taken with a different realization of the trap de- 



scribed in Sec. |IV A| but of the same design, so imperfec- 
tions in construction and small variations in the electrode 
impedance may have affected the measurement. Differ- 
ences in impedance affect the strength of the gradient, so 



the agreement in B' values is probably fortuitous. How- 
ever, we find again a value for aM W that is close to the 

-25° (Sec.|lVA[). 



simulated value 



C. Control of the orientation of the radial motional 
modes 

After determining the characteristics of the microwave 
quadrupole, we align the orientation of the radial mo- 
tional modes to the microwave field gradient to op- 
timize sideband Rabi rates. Since the qubit is con- 
trolled with a Amp = transition, its motional side- 
bands are driven by 5B\\/8r (Fig. [8^1). To overlap 
the HF-mode vector with 5B\\/5r\ max , we use an ad- 
ditional control potential (j) TOt with potentials Vq\_ cq = 
{-1.694,2.298,-1.694,0.311,0.501,0.311} V applied to 
the control electrodes, which leads to a mode orienta- 
tion «lf — 5° in the absence of any additional con- 
trol potential. In the experiments, the total potential 
0totai = </>tra P + Astray + s TOt (j) TOt breaks the symmetry of 
the rf pseudopotential, and sets the orientation of the 
radial modes. The value of the scaling coefficient s ro t is 
used to adjust the orientation of the radial modes. For 
the given total potential </>totai including the stray po- 
tential (Appendix [A]), negative ,s ro t rotates the LF-mode 
vector toward «lf —> 0° and for positive values it rotates 
toward q^lf —> —90°. Furthermore, rot is designed such 
that the ion position remains unchanged while all trap 
frequencies change as a function of the scaling. A com- 
parison of the simulated frequency changes with the ac- 
tual experimental data can be used to determine and ad- 



just the orientation of the radial modes. In Fig. 10 1 such 
a comparison between simulated and measured mode fre- 
quencies shows good agreement. Figure [Tofc illustrates 
the corresponding mode-angle values extracted from the 
simulation. Note that for the simulation we had to in- 
clude Astray to find agreement with the experimental 
data, details are given in Appendix |Aj When choosing 
«lf = 2ctMW — —53°, we observe the maximal cou- 
pling of the microwave-field gradient to the HF mode, 
while the LF mode couples very weakly to the microwave 
field. We find a corresponding Rabi rate for driving the 
HF sideband transition of the motional ground state of 
~ 2ir x 2.0 kHz, with the microwave-field gradient ad- 
justed to 35 T/m. 



V. PRESENT LIMITATIONS FOR TWO-QUBIT 
GATES 



Internal-to-motional state coupling as well as an en- 
tangling two-qubit gate based on the simultaneous ap- 
plication of detuned red and blue sideband drives [42] 
was presented in [32]. Using this gate, entangled states 
with a fidelity of approximately 0.76 were produced. 
In these experiments we used an out-of-phase radial 
mode of two ions, because it is less sensitive to electric- 
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FIG. 10. Control over the radial mode orientation, (a) 
Comparison of simulation (lines) with experimental data of 
radial-mode frequencies (error bars are too small to be vis- 
ible) as a function of the scaling of the control field, (b) 
Mode angle «lf extracted from the simulation. At an an- 
gle «lf = 2qmw — —53° the high-frequency mode is aligned 
along 5B\\/ 5r\max and couples maximally to the microwave- 
field gradient. 



field noise than the corresponding center-of-mass mode. 
The out-of-phase motional heating was measured to be 
0.2 — 0.5 quanta/ms, while the center-of-mass motional 
heating rate was approximately 5 quanta/ms (for radial 
mode frequencies around 7 MHz). Heating correspond- 
ing to an increase of a single quantum of motion on the 
mode used for an entangling-gate operation effectively 
decoheres the gate operation [52 j. Since in [32] the en- 
tangling gate durations were approximately 250 /is, even 
the out-of-phase mode heating can cause substantial gate 
errors. To investigate this effect for our experimental con- 
ditions, we numerically solve the master equation given 
in [52]. We determine that this effect alone limits the 
fidelity to approximately 0.8. Furthermore, we observed 
drifts of the motional sideband frequencies in the exper- 
iments, which may be caused by mode frequency drifts 
and/or varying ac Zeeman shifts from drifts of the resid- 
ual magnetic and electric microwave fields. We estimate 
the frequency variation during the gate pulse to be less 
than 0.5 kHz from these effects, while a typical detuning 
for the gate is around 4.5 kHz. This systematic effect can 
contribute an additional 5 % loss of fidelity. 



VI. CONCLUSION 

In summary, we have given details of the surface- 
electrode trap and microwave electronics used for the 
microwave near-field quantum control experiments de- 
scribed in [32]. We have described the experimental 
procedures used to control, shape and calibrate the mi- 
crowave fields that implement motional-sideband tran- 
sitions using the near-field gradients. We have intro- 
duced a novel technique for rf micromotion detection. 
Furthermore, we presented and experimentally verified a 
method to rotate the radial mode axes. Finally, we in- 
vestigated limitations in fidelity of the two-qubit gates 
in [32] through simulations, and examined the impact of 
leading error sources. We determined that the motional 
heating during the gate pulse is the dominant source of 
infidelity. A significant reduction (by a factor of 100 or 
more) of the heating rate should lead to entangled state 
fidelities higher than 0.9. A reduction of ambient heating 
may be achieved by surface treatments of the electrode 
structure as described in [38[|39] and/or cryogenic cooling 
of the electrodes [35-37 . In addition, higher microwave 
currents can increase the gate speed and thereby can re- 
duce the effect of motional heating. The corresponding 
increased heat load may be mitigated by reducing elec- 
trode resistance through low temperature operation. 
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Appendix A: Comparison of Trap Simulation with 
Measurements 



To compare the trapping potential calculations from 
Sec. [H] to the experiment, we first determine the electric 
field strengths at the p ositi on of the ion to compensate 
the rf micromotion (Sec. IV). We then measure the mode 
frequencies, and find the orientation of the radial mode 
vectors (Sec. IV C), characterized by «lf (Fig. [2]). We 
find a stray field E stray ~ 720 x {-0.54, 0.10, -0.84} V/m 
and a stray field curvature tensor, given by the Hessian 
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matrix of (/) s t ray at the center of the trap, 

1.96 0.44 0.39 
Astray ^ ( 0.44 -1.67 -4.23 ) x 10 7 V/m 2 
0.39 -4.23 -0.29 

After including </> s tray into the potential simulations, we 
find reasonably good agreement (motional frequencies 
and mode angles agree to within 10 % or bette r) be- 



tween the simulations and the measurements (Sec. IV C ) 
for several different settings of rf voltage and control po- 
tentials. The stray potential might be caused by work 
function variations on the electrode surface [53H56] (e.g., 
from varying Mg contamination) or stray charges local- 
ized in the central region of the trap chip. One way of 
modeling modeling the stray effects is a rectangular patch 
with size Ay = 35 fim and Az = 5 /im located on the 
chip surface, whose center is displaced by z = —25 /im, 
y = +3 /im relative to the trap center and that is biased 
to —1.15 V (typical work function differences can exceed 
one volt). We observe temporal variations in the stray 
potential effects, perhaps due to uncontrolled charging 
of the electrodes which may be caused by laser light im- 
pinging on the electrodes and/or changing Mg coverage 
of the chip surface during loading of the trap. In partic- 
ular, the applied compensation fields vary within a few 
tens of volts per meter on a day-to-day basis, and the 
axial-mode frequency can vary more than 200 kHz over 
the course of a few weeks with many loading cycles. 

We estimate the deviation of ion displacements from 
simulated values due to stray-electric fields to be less than 
a few percent for typical displacements from the rf null of 
±1 /im at radial- mode frequencies of ~ 7.0 MHz. We do 
not include this effect in our analysis of the experiments 
described in Sec. |IV A| and Sec. IV B| However, for the 
calibrations described in Sec. IV C[ we take the stray- 



electric field into account because of its significant impact 
on the radial motional-mode vector directions. 



Appendix B: Influence of the microwave electric 
fields on the trapping potential 

In addition to the oscillating magnetic field Bmw, an 
electric field Emw ( a t the microwave frequency) is also 
present (Fig [8)3). If the null points for Emw and Bmw 
fields differ, the pseudopotential produced by Emw can 
push the ions away from the main trap rf pseudopo- 



tential null. From numerical simulations, we determine 
that the ion position in a 7.0 MHz trap is shifted by 
{5x, Sz} ~ {0.3, -0.2} nm when B f ~ 35 T/m (Fig. |]o). 
This displacement is much smaller than the precision for 
positioning the ion with the shim fields (Sec. [H| and 
therefore this shift cannot be resolved. 

Another systematic effect is caused by the curvature 
of the microwave pseudopotential. It modifies the mode 
frequency (based on simulations) by about 0.5 kHz, de- 
pending on the mode direction. This can be a significant 
shift for the two-qubit gate. However, when calibrat- 
ing the sideband frequencies this effect is included, since 
we probe the sidebands with the microwave drive, which 
inherently includes the microwave pseudopotential pro- 
duced by the gate pulse. To experimentally check for 
this shift, we measure the mode frequencies by applying 
a oscillating potential (coherent excitation) to one of the 
control electrodes, while the ion is in the |3, 3) state. We 
then detect the ion fluorescence as a function of excita- 
tion frequency: a decrease in the resulting fluorescence 
indicates that a mode of the ion's motion has been res- 
onantly excited [57]. We compare the measured mode 
frequencies with and without the microwave fields while 
applying the excitation pulse. With the microwaves ap- 
plied, we observe shifts of the radial- mode frequencies of 
~ 3 kHz. The observed shift is about seven times higher 
than what we expect from our estimate; this discrep- 
ancy is not understood. The mode frequency shift (for a 
7.0 MHz trap) for the ion in the |3,3) state due to the 
spatially varying ac Zeeman shift is estimated to be less 
than ~ 10 mHz and can be neglected here. 

We checked for motional excitations when simultane- 
ously applying two tones near /o ± /radial > as required for 
the two-qubit gate (Sec.[V]). For pulse lengths up to 1 ms 
(four times longer than typical gate times) we observe no 
change in the motional state within an uncertainty of ap- 
proximately ±0.5 quanta. However, when we apply the 
two tones at fo ± / ra diai/2, we record coherent motional 
excitations corresponding to a few tens of quanta after 
50 /is. The origin of the observed excitation may be the 
gradients of the two microwave pseudopotentials (with 
difference frequency /radial)- We estimate this effect to 
be negligible for our experimental conditions. However, 
we cannot rule out that the excitation is caused directly 
by an electric field at frequency /radial- This field may 
originate from mixing of the two microwave tones due 
to nonlinear behavior of one or more of the microwave 
elements. 
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